Identified Particle Spectra for Au+Au Collisions at y/s = 200 GeV from 
STAR, PHENIX and BRAHMS in Comparison to Core-Corona Model Predictions 
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The core-corona model describes quite successfully the centrality dependence of multiplicity and 
< pt > of identified particles observed in heavy ion reaction at beam energies between ^/s = 17 
GeV and 200 GeV. Also the centrality dependence of the elliptic flow, V2, for all charged and 
identified particles could be explained in this model. Here we extend this analysis and study the 
centrality dependence of single particle spectra of identified particles. We concentrate here on 
protons, antiprotons, kaons and pious which have all been measured by the PHENIX, STAR and 
BRAHMS collaborations. We find that an analysis of the spectra in the core-corona model suffers 
from differences in the data published by the different experimental groups, notably for the pp 
collisions. For each experience the data agree well with the prediction of the core-corona model but 
the value of the two necessary parameters depends onthe experiments. 



I. MOTIVATIONS 



Lattice calculations predict that at high density and temperature hadronic matter transforms into a plasma of 
quarks and gluons. There is evidence that such a state can be obtained in heavy ions collisions at beam energies 
which can be reached at the colliders at CERN and in Brookhaven. This plasma is a very short living state - it lasts 
^ ] less than 10^^"^ seconds - but it is assumed that this time is sufficiently long for reaching equilibrium. Hydrodynamical 
calculations indeed shovif that the experimental spectra can be described assuming that an equilibrium is obtained in 
a very short time. The origin of such a fast equilibrium is still discussed. The matter expands quickly toward the 
chiral/confinement phase transition during which hadrons are formed. It came as a surprise that the multiplicity of 
identified stable particles in the most central collisions agrees almost perfectly with that expected for a statistical 
distribution at a freeze out temperature of around 170 MeV and a small baryon chemical potential. This has been 
taken as well as evidence that at least at the chiral/confinement phase transition the system is in equilibrium. 
] For symmetric systems the number of projectile participants equals that of target participants, independent of 
fS| ' the centrality. If each particpant contributes the same energy in the center of mass system and if the system come 
I to equilibrium one does not expect that the multiplicity per participating nucleon varies with centrality. In the 
I i experiments such a variation has been observed, however. In addition it depends strongly on the particle species. 
' Whereas for tt this ratio is almost constant, for multi strange baryons this ratio varies by a large factor, a phenomenon 
^-H . which has been dubbed strangeness enhancement. 

^ ' One of the crucial assumptions in statistical model calculations is that geometry docs not play a role and all nucleons 
. ^ , come to statistical equilibrium. In simulations of the heavy ion reactions it has been observed that especially those 
' nucleons which are at the edge of the overlap region have initially few, in some cases only one, collision [l| . It does not 
H I seem to be realistic that those equilibrate. The ratio of these edge nucleons as compared to all participating nucleons 
- - 1 decrease with centrality. They do not present an important fraction for the most central collisions. This observation 
has motivated the core-corona model in which is it assumed that nucleons which scatter initially only once (corona 
particles) are not part of the equilibrated source but produce particles as in pp collisions whereas all the other come 
to statistical equilibrium (core particles) . Of course this fast transition between core and corona particles is a crude 
approximation but it allows to define from experimental pp and central AA data the centrality dependence of the 
different observables. Studies have shown that the present quality of data does not allow for a more refined definition 
of the transition between core and corona particles. It has been further verified that the core-corona model describes 
quantitatively the results of the much more involved EPOS simulation program. 

In a series of papers [H-Q it has been shown that the core-corona model describes quite nicely the centrality 
dependence of the multiplicity of identified particles, < pt > of identified particles and even of V2 observed in AuAu 
and PbPb collisions. The latter has been considered as a test ground for the shear viscosity needed to describe heavy 
ion data in viscous hydrodynamical calculations. The core-corona model describes this data without any reference to 
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a viscosity. The prediction for the CuCu data are completely determined by the AuAu data and agree with data as 
far as data have been published. 

In this contribution we go one step further and study the observed single particle pt spectra for antiprotons, protons, 
kaons and pious. This study has the following objectives: 

a) To investigate whether the core-corona model, which has been successfully applied to described the centrality 
dependence of several averaged variables of identified and not identified particles like the multiplicity, < pt > and the 
elliptic flow V2 , is also able to describe the centrality dependence of the spectra 

b) Experimentally is has been observed that the centrality dependence of < pt > for tt, K, p and p is rather different. 
This has been interpreted as a sign for collective radial flow which yields a stronger increase of < pt > with centrality 
for heavy hadrons as compared to light ones. Here we address the question whether this increase is really only due to 
the superposition of core and corona spectrum or whether the increase of the mean values signals deviations from a 
simple superposition. 

c) To search for domains in pt which show deviations from the core-corona prediction, expected if interactions between 
core and corona particles become important, and to try to interpret those deviations, if they exist, in physical terms. 

II. CENTRALITY AND CORE-CORONA FRACTION 

In order to determine the centrality dependence of the specta we have first of all to know the relative contribution 
of core and corona particles as a function of the centrality. The core-corona model relies on a single parameter : 
f{Ncore), the fraction of core nuclcons as a function of the centrality. Along with the number of participans, Npart, 
it is calculated by a Monte-Carlo simulation based on a Glauber model for hadrons in the nucleus. The parameters 
of the Glauber distribution are the only freedom of this model. We apply here the EPOS approach. The results are 
presented in table 1. We present there for the different centrality bins the average number of participants. Npart, as 



Centrality 


80%-92% 


70%-80% 


60%-70% 


50%-60% 


40%-50% 


30%-40% 


20%-30% 


10%-20% 


5%-10% 


0%-5% 


^ ^ part 


7.939 


16.40 


29.87 


51.01 


80.32 


119.0 


169.1 


235.1 


300.9 


351.3 


fcore ^Npart ) 


0.269 


0.407 


0.504 


0.587 


0.669 


0.727 


0.782 


0.831 


0.869 


0.911 


mSTAR 
''-^ part 




10.7 


22.0 


40.6 


67.8 


105.4 


155.9 


223.6 


289.6 


345.8 


a.^tSTAR 




2.1 


3.2 


4.3 


5.0 


5.3 


5.1 


4.2 


3.0 


1.9 


jyrPHENIX 
^' part 


6.3 


13.4 


25.7 


45.5 


74.4 


114.2 


166.6 


234.6 


299.0 


351.4 


A ktPHENIX 


1.2 


3.0 


3.8 


3.3 


3.8 


4.4 


5.4 


4.7 


3.8 


2.9 



TABLE I: Centrality dependence of f{Ncore) and Npart for STAR, PHENIX and the core-corona model. 



well as the core fraction fcore- The STAR g] and PHENIX Q collaboration found other and mutually different values 
of Npart for the same centrality bins. The values of the two collaborations agree within error bars but since this is 
a completely theoretical quantity it is not clear why the difference cannot be avoided. Because different Npart yield 
different fcore it becomes difficult to compare the different experiments with the same model parameters. 

III. RESULTS 
A. Multiplicity 

In this contribution we investigate the centrality dependence of the tt, K and proton data measured by all three 
RHIC collaborations, STAR Q, PHENIX Q and Brahms @,[l3l • According to the publications statistical errors are 
small and therefore systematical errors dominate. The spectra of tt and K of the three experiments should be directly 
comparable. For the proton spectra two collaborations (BRAHMS and PHENIX) have corrected the spectra for feed 
down due to weak decay whereas STAR published uncorrected data. This is regrettable in view of the fact that the 
vertex resolution of STAR is better than that of PHENIX. Brahms, however, did not apply a feed down correction 
for the preliminary proton spectrum in pp collision yet [loj . In order to apply the core-corona model to the Brahms 
data we have assumed that 40% of the protons in pp are due to weak decay, employing the reduction factor published 
by the PHENIX collaboration [|. 

Fig. 1 displays the experimental data. On the left hand side we have divided the published value of the multiplicity 
by Npart calculated in the EPOS approach for the concerned centrality bin (see table 1). On the right hand side we 
have divided the multiplicity by the Npart value which has been published by the experimental groups. We observe 



3 




• 


P 


■ 


p (xO.3) 


♦ 


A'+ (x2) 


A 




9 


STAR 


• 


BRAHMS 


• 


PHENIX 


• 


Data 




— Core corona 



^ 10-' li 





e^s^ ^ ^ 

mo OOOOOOO tN 

^^^^^^^ ^ 

mOOOOOOO O 



(a) Core - Corona Npart 



(b)Respective experimental Npart 



FIG. 1: Multiplicity per participant as a function of centrality. On the left hand side we use Npart as calculated by EPOS, on 
the right hand side the different Npart quoted by the experiments for the given centrality bin. Symbols are data, the core-corona 
predictions are displayed as lines. 



that for TT and K the multiplicities for central events agree quite well and are in the mutual error bars but if one looks 
closer one realizes that STAR observes a increase of 10% of the multiplicity of the pious toward the most central 
bin which is not seen by PHENIX. Toward peripheral collisions the differences increase and the data are outside 
the mutual error bars. For the multiplicities observed for pp resp. very peripheral AA {Npart=Q.9) collisions we 
observe differences of 30%. The STAR data set is almost compatible with a centrality independent tt multiplicity per 
participant, whereas for the PHENIX data set this is certainly not the case. If one takes the different Npart values of 
the different experiments, right, the differences become smaller. 

For p and p the heavy ion data of BRAHMS and PHENIX agree but the preliminary pp data points of BRAHMS 
are almost a factor of 2 higher than that for the most peripheral data point of PHENIX. In ref. Q the PHENIX 
collaboration quotes that between 30% and 40% of the observed protons and anti-protons are due to weak decay. 
Therefore it is surprising that the multiplicity observed by STAR and PHENIX differ by more than a factor of 2 
in central collisions and by an even larger factor (close to 3) for peripheral reactions where the correction due to 
weak decay should be smaller (strangeness enhancement). Also the trends are quite different. Whereas STAR and 
BRAHMS observe a quite modest centrality dependence PHENIX observes a quite strong one. 

In the core-corona model the centrality dependence of the multiplicity of a given particle species i in a centrality 
bin containing Npart participants is given by : 

M'iNpart) = Npart. [fcore{Npart) -M'^ore + (1 " fcore{Npart)) -M'^orona] (1) 

where M^^^^ is the multiplicity per core participant and M^^^^^a the multiplicity per corona participant. There are 
several ways to determinate these two values: one can either calculate them from integrated fits or one can use directly 
the published values, which are the results from a fit of a specific form (blast wave model) to the experimental spectra. 
We chose the latter in the present study: we use the multiplicity measured in pp and divided by a factor of two for 
-^corona- Then wc extract M^^^,^ from the most central multiplicity using eq. (1). If there was no suitable pp data like 
for PHENIX, we use instead the most peripheral AA bin and eq. (3) to determine Ml^^^^a- The result is displayed 
in fig. 1 as straight line. We observe that for all experiments the centrality dependence is that which is expected in 
the core-corona model but, as already said, the large difference between the experiments does not allow for an unique 
value of M*„^„„„ and of M^^^^. 



B. Spectra 



It is generally believed that the hadrons, after being formed during the confinement phase transition, interact on 
the way to the detector. Event generator are based on this assumption but the results of EPOS [1] demonstrates that 
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this rescattering changes the spectra only at low pt, where almost no experimental data are available. EPOS succeeds 
to reproduce the measured spectra in between a factor of two and, in consequence, to reproduce the change of the 
spectral form from central to peripheral collisions. The origin of this success, however, is little transparent due to the 
complexity of the approach. Therefore it is worthwhile to study the spectra in an approach which contains the (for 
the success of EPOS) essential separation of core and corona particle but which is otherwise as simple as possible, 
i.e. in the core-corona model [2tJ1- this model we can calculate which spectra would be expected if no final 
state interactions among hadrons take place and we can use the difference to the data to learn something about the 
final state interaction. Assuming no final state interactions, in the core-corona model the spectra are superpositions 

^corona 

of two contributions: the core contribution and the corona contribution. The corona distribution 7^ S — j- is the 

^ 2-ir.pt.dpt.dy 

measured pp spectra divided by two, the core contribution 2tt p^dpt dy ^® obtained from the experimental spectra for 
the most central AA collisions corrected for the corona contribution (sec tabic I) and divided by the number of core 
participants. Then in the core-corona model the spectra for a given centrality is given by: 

cPA'I- ^2 p^corona ^2 j\^core 
7; 1 -f- = Npartii'^ ~ !core{J^part))-Z 1" + fcore{Npart) 7: ^ 1"]- (2) 

2Tr.pt.dpt.dy 2Tr. pt.dpt.dy 2t:. pt.dpt.dy 

Fig. [2] presents the results for mesons. We see in the top row the spectra measured by the different experimental 
groups in comparison with the predictions of the core corona model. The left hand side presents the STAR data, the 
middle panel the PHENIX data and the right hand side the BRAHMS data. The middle row displays the difference 
between theoretical predictions and the data. The shaded region mark the error bars (which are taken as the averaged 
error bar over the centrality bins). We observe that for almost all STAR data points the theoretical predictions are 
in the experimental error bars. In peripheral reactions there is a tendency that at large pt the core-corona model is 
above the data. For semi central reactions model and data are in agreement for all pt values. This is not trivial at 
all as the bottom row shows. There we display the inverse slope parameter obtained by fitting the experimental and 
theoretical spectra by a thermal spectra 

— = C.mt.e— 3 

Zir.mt.dmt.dy 

with rrit being the transverse mass of the considered particle: 
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y/jH^+m? (4) 



where pt is the transverse momentum of the particle (with respect to the beam axis) and m its free mass. Even if the 
curves are not exactly exponential and therefore the value of the inverse slope parameters depends on fit range the 
bottom panel shows clearly that the slope varies considerably from central to peripheral reactions (which is in the core 
corona approach a consequence of the different invariant slope parameters in pp and central AA collisions) . Also the 
PHENIX data are compatible with the core-corona model besides the second last and third last centrality bin where 
the deviations, expected from fig. 1, show up. We note in passing that for those centrality bins the core-corona models 
agrees well with the STAR data. The central BRAHMS data are also compatible with the model but we see deviations 
for the most peripheral bin. It is remarkable that the peripheral STAR data are almost exponential whereas those of 
PHENIX arc not and even less those of the BRAHMS collaboration. Comparing the three experimental spectra with 
the model wc can conclude that for each experiment the majority of data is well described by the model. Deviations 
are specific for the experiment. There are no systematic deviations. 

Fig. [3] show the same quantities for the protons. For the STAR data see an almost perfect agreement between 
data and model predictions. The only exceptions are, as for the K~^, peripheral data at large pt- The middle left 
figure demonstrates this in detail. We see that besides the high pt points in peripheral collisions the spectral form 
is reproduced by the core-corona model in between the error bars. This is far from being trivial. The inverse slope 
parameter of the spectrum varies by a factor of 2 between central and peripheral reactions and so does the inverse 
slope parameter of the model due to the large difference between in the inverse slope parameters in pp and central AA. 
Almost the same is true for the PHENIX data. Here the model ovcrpredicts the data of the second to last centrality 
bin by an almost constant factor. That this bin is particular we have already discussed in fig.l. In the BRAHMS 
proton data the form varies from central to peripheral reaction, what is not seen in the STAR and PHENIX data. 
This can also not be reproduced in the core-corona model. We do not display here the comparison for the tt, which is 
also reproduced in the core-corona model, because neither the multiplicity per participant (see fig. 1) nor the slope 
varies considerable from central to peripheral reactions. Therefore it is not a good test ground for the validity of the 
core-corona model. 

Thus the and p spectra can be described in the core-corona model, however with rather different distributions 
for the core and corona particles. Why different parameters are necessary is demonstrated in detail in figs. U and [5] 
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FIG. 2: Centrality dependence of the K+ spectra measured by the STAR (left), PHENIX (middle) and BRAHMS (right) 
collaborations in comparison with the predictions of the core-corona model. Top: The experimental spectra in comparison with 
the model prediciton. Middle: Ratio of the prediction of the core-corona model and the experimental data. The shaded area 
mark the experimental errors. Bottom: Inverse slope parameter T obtained by fitting data and theory by eq. 3. 



in which wc compare directly the data of the PHENIX and of the STAR collaboration for a fixed value of pf . Fig|3] 
shows how the yield of this pt bin develops as a function of the centrality. Fig. [S] displays the ratio of the yields 
measured by STAR and PHENIX as a function of the centrality. 

For the tt"*" and the the form of the spectra measured by STAR and by PHENIX deviate for peripheral reactions 
and for the in addition also for the most central collisions. Therefore the ratio of the spectra deviates from the 
expected value of 1 ( shown as a straight line in fig. [5]). Only three out of the 18 points agree in the error bars ( 
which are added quadratically) with 1. 

The proton and antiproton spectra should not be identical because the PHENIX data have been corrected for weak 
decay whereas the STAR data are not. According to the PHENIX results this correction should be around 30 % 
- 40% . This would explain a ratio of 1.66 between the data shown as well as a straight line in fig [SJ The data 
differ, however, by more than a factor of two for central and semi central reaction and the difference increases toward 
peripheral reactions, as seen in fig. [5] 
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FIG. 3: Centrality dependence of the proton spectra measured by the STAR (left), PHENIX (middle) and BRAHMS (right) 
collaborations in comparison with the predictions of the core-corona model. Top: The experimental spectra in comparison with 
the model prediciton. Middle: Ratio of the prediction of the core-corona model and the experimental data. The shaded area 
mark the experimental errors. Bottom: Inverse slope parameter T obtained by fitting data and theory by eq. 3. 



In conclusion we have shown that all three available data sets for pt spectra of protons and kaons can be well 
described in the core-corona model. There are deviation but where they occur varies from experiment to experiment. 
Unfortunately no common parameter set can be found which describes the three experiments simultaneously. This is 
due to the strong differences between the experimental results for more peripheral events and due to the differences 
of the spectra measured in pp where the multiplicities differ up to a factor of two between the three experiments. 

The fact that the core-corona model describes also the single particle spectra of identified particles has important 
consequences: 

a) The fact that the spectra at different centralities are only a superposition of that of core particles which are in 
thermal equilibrium and that of corona particles which have a spectrum as measured in pp collision testifies that 
there is little interaction between core and corona particles. Otherwise momentum would be transfered between the 
different particle species which have a quite different average pt- This would, in consequence, modify the spectral 
form. 
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FIG. 4: Comparison of of the centrality dependence of the yield measured by STAR and PHENIX for a given pt- 
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FIG. 5: Ratio of the yield for a given pt measured by STAR and PHENIX as a function of the centrality. 



b) The observed different centrality dependence of the average pt for different particles which has been interpreted as 
a signal for collective flow is merely due to the fact that the difference of pt in central AA and pp collisions is particle 
specific. The difference is large for protons, practically inexistent for pions and that for kaons is in between. Of course 
it may be that the collective flow is the origin of the different values for pt in central heavy ion and pp collisions. The 
core-corona model tell us that then the collective flow is independent of the size of the core part. 

c) The blast wave fit, usually applied to characterize the data, is not the proper tool to analyze the spectra. Rather 
then describing a two component spectra by one single fit it would certainly be worthwhile to subtract the corona 
part in oder to study the properties of the core part which may carry the desired signal of the plasma. 

Adding these observations to the already published analysis of the multiplicity HI , < > and V2 of identified particles 
[3, 01 in the core corona model one can conclude that presently all data on the centrality dependence of observables 
can be quantitatively described in this model. Predictions for the CuCu are straight forward and do not need any 
further input. It remains to be seen whether other approaches based on non-viscous hydrodynamics yield a similar 
good description of the data. This would a prerequisite for validating the idea that the centrality depedence of W2 /e 
is really due to a finite viscosity and not due to a distinction between core and corona particles. 
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and for an interesting discussion. 
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